Finite-rate quenches of site bias in the Bose-Hubbard dimer 
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For a Bose-Hubbard dimer, we study quenches of the site energy imbalance, taking a highly 
asymmetric Hamiltonian to a fully symmetric one. The ramp is carried out over a finite time that 
interpolates between the instantaneous and adiabatic limits. We provide results for the excess energy 
of the final state compared to the ground state energy of the final Hamiltonian, as a function of 
the quench rate. This excess energy serves as the analog of the defect density that is considered in 
the Kibble-Zurek picture of ramps across phase transitions. We also examine the fate of quantum 
'self-trapping' when the ramp is not instantaneous. 



I. INTRODUCTION 

Explicit time evolution of quantum many-particle sys- 
tems out of equilibrium has generated intense interest 
due to unprecedented possibilities for experimentally ac- 
cessible non-equilibrium situations, opened up by devel- 
opments in laser-cooled atomic clouds and in mesoscopic 
systems. One theme has been the response to an in- 
stantaneous 'quench', where a physical parameter is sud- 
denly changed to a different value. While instantaneous 
quenches are more convenient to analyze, a change of pa- 
rameter can of course be performed at any rate. At the 
other extreme from the instantaneous quench, one can 
make the parameter sweep adiabatically, in which case 
the system reaches the ground state of the final Hamil- 
tonian. The Kibble-Zurek theory describes finite-rate 
ramps, neither instantaneous nor adiabatic, across phase 
transitions ji While the original interest concerned ther- 
mal phase transitions, there has been a recent surge of 
interest in finite-rate traversals of quantum phase tran- 
sitions in lattice Hamiltonians (Refs. and references 
therein). One important context is the Bose-Hubbard 
Hamiltonian, where finite-rate ramps of the Hubbard in- 
teraction U have been analyzed^ - — motivated by an influ- 
ential experiment quenching across the Mott-supcrfluid 
transition^ 

In this work, we analyze finite-rate ramps in a Bose- 
Hubbard dimer. For far-from-equilibrium issues where 
few standard theoretical techniques exist, it is of obvi- 
ous interest to look at finite clusters because of near- 
exact solvability. One can hope to understand the non- 
equilibrium dynamics in some detail. Such detailed 
results for clusters clearly provide an invaluable back- 
ground for the emerging field of non-equilibrium dynam- 
ics in macroscopic (many-site) quantum systems. In ad- 
dition, the bosonic dimer is an important model sys- 
tem by itself, and its physics is relevant in diverse con- 
texts. In recent years several experiments have achieved 
two-site bosonic systems in cold-atom setups^ - — The 
quenches we study could be implemented in such a 
setup, e.g., through slower ramps in the experiment of 
Ref. H. The Bose-Hubbard dimer can be mapped onto a 
single-spin Hamiltonianji 3 - - — very similar to that govern- 
ing single-molecule magnetic experiments, where finite- 



rate ramps have received some attention! 20 ' 21 Finite-rate 
quenches have also been reported with a Josephson junc- 
tion arrangement^ 2 - which is closely related to a bosonic 
dimer. 

For N bosons in two sites, the Hamiltonian is 
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with ni = a\ai and (raj.) + (n 2 ) = N. 

Quenches of the interaction U in the dimer have been 
considered in Refs. 4-6.23. However, for the dimer, 
quenches of the site bias parameter 8 are also of natu- 
ral interest. Instantaneous quenches of 8 from an un- 
balanced to a symmetric situation, i.e., S(t) = 8q6{—£), 
have been studied as a quantum realization of self- 



trapping. 
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Classical self-trapping, described by the 



two-mode Gross-Pitaevskii equation, is the persistence of 
imbalance in the state despite the final Hamiltonian be- 
ing unbiased. 8 'i21 In the quantum case, the relative num- 
ber imbalance z = (nj — n2)/N oscillates after the in- 
stantaneous quench, with decay and long-time revivals 
of the oscillations. For small So, the imbalance oscilla- 
tions are around zero. However, for large enough values 
of Sq and NU/K, after a fast quench the mean value 
of the oscillating z(t) does not relax to zero within any 
reasonable time scaleJ^r— The astronomically long re- 
laxation times ("quantum self-trapping") can be under- 
stood in terms of extremely small energy splittings in 
the eigenvalue spectrum of the final Hamiltonian. For 
the instantaneous quench, the various relationships be- 
tween timescales of dynamical features and energy scales 
of the 8 — spectrum has been discussed in some detail 
in the literature.— ~— 

We use U = 1, measuring energy [time] in units of U 
[h/U]. We focus mostly on the regime K «U. We will 
analyze finite-rate quenches of 8, of the form 

8{t) = -So 6{-t) + (-8oM(T-t) 9(t) Q(t — t) , 

as illustrated in Figs. UK, [lb. We consider the whole 
range from r = (instantaneous) to r — > oo (adiabatic). 
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FIG. 1: (Color online.) (a) and (b) Form of quench ana- 
lyzed in this work. Lower panels illustrate the approach to 
adiabaticity with increasing r, for (N,K,5o)=(i, 0.2, 9). (c) 
overlap with ground state of final Hamiltonian. Curves from 
lowest to highest final values correspond to r = 10, 40, 80, 
150. (d) Evolution of energy: curves from top to bottom are 
r = 0.01, 10, 40, 80. Dotted horizontal line is the ground 
state energy of the final Hamiltonian. The final excess energy 
over this dotted line is the residual energy, which decreases 
with t (e). 



The initial (t = 0) state is taken to be the ground state of 
the initial Hamiltonian with 5 = —So- The initial asym- 
metry So is taken to be large enough that the bosons are 
initially concentrated almost entirely on site 1, i.e., the 
wavefunction is dominated by | iV, ) . The final ground 
state is dominated by \N/2,N/2) but the system does 
not reach this ground state unless the quench is truly 
adiabatic. Fig. [TJ; demonstrates how larger-r quenches 
are more nearly adiabatic, through the temporal evolu- 
tion of overlaps with the 5 = ground state. 

In Kibble-Zurek theory, one uses the defect density in 
the final ordered state to quantify the deviation from 
adiabaticity. In few-site clusters, defects are not easily 
defined, nor are phase transitions or ordering. How- 
ever, there is a natural quantity that serves an anal- 
ogous role, namely, the final energy after the quench, 
Er>T = {ipt>r \ H s=0 |^>t> T ). In an adiabatic sweep, the 



final energy is the ground state energy E^ of the fi- 
nal Hamiltonian. The excess energy over Eg^ (residual 
energy) measures the deviation from adiabaticity - The 
interpolation between instantaneous and adiabatic limits 
is illustrated through energy evolution in Fig. [TJl, and 
through residual energies in Fig. [T£. 

Our main results concern the dependence of the resid- 
ual energy, AE = E t>T — E s g = & ° , on the quench time r. 
For near- instantaneous quenches (small r), the energy 
deviation from the instantaneous limit is found to scale 
as t 2 . We analyze larger-r quenches through a multi- 
crossing Landau- Zener— scenario, and derive concise ex- 
pressions for AE in an intcrmediate-r regime as well as 
in the near-adiabatic (very large r) limit. 

While we focus on AE, our analysis can in principle 
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FIG. 2: (Color online.) (a-c) energy spectra for N = 6 
bosons. For small K and even N, the lowest state goes 
through N/2 avoided level crossings on either side of 8 — 0, 
at S = ±{2x - 1)U, with x = 1, 2, N/2, (d) N = 2 bosons, 
(e) Residual energy; dots are exact values and line is theory. 



be adapted to treat other observables. Since instanta- 
neous (5-quenches are associated with the self-trapping 
phenomenon, we also ask whether and how self-trapping 
survives when the ramp rate is finite. 



II. SLOW QUENCHES 

Our analysis for large r relies on the avoided level 
crossing structure of the problem. Fig. [5^-c shows that 
the level crossing structure is only relevant for small K. 
Our analytic treatment for slow quenches is focused on 
this parameter region, K <C U ( Fig. f2fi). Therefore, 
we can regard the quench problem as one of traversing a 
sequence of well-separated avoided level crossings. This 
scenario would clearly not be applicable for situations 
such as those in Figs. [2b, Ht. 

We start with the simplest case of two bosons. There 
is only one level crossing encountered during the quench, 
at 5 w — 1, where the |2,0) and |1,1) states are mixed 
(Fig. The Hamiltonian within this space is H = 

1 + S(t) K/^/2\ ^ Landau-Zener formula 25 gives 
-K/V2 J 5 
the probability of excitation at this level crossing to be 
p = er 2 ^, where 7 = (K/V2) 2 /S. At the end of the 
quench, the energy has a contribution of weight (1 — p) 
from the ground state (energy sa 0) and a contribution 
of weight p from the two higher levels (energy sa U). 
Thus the residual energy is AE sa Ue^^ /So)T . This 
reproduces numerical calculations for large r (Fig. [2^) . 

For larger numbers of bosons, the quantum state can 
take various paths to 6 = 0, in a multi-crossing situa- 
tion such as that shown in Fig. [2^,. Fortunately, at small 
K the higher crossings are simple to treat because the 
energy splitting at these points are of higher than lin- 
ear order in K. Thus for K<tiU one can regard these 
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as real crossings rather than avoided crossings, so that 
the excitation probabilities are unity. One therefore has 
to consider only the excitation probabilities at the N/2 
crossings involving the lowest energy state. The excita- 
tion probabilities are 

Pi = exp[-f i{N - i + 1)K 2 t/5 ] 

at the i-th crossing encountered during the quench; here 
i(N — i + 1) is the bosonic factor relevant for the coupling 
between states \N — i + 1) and \N — i,i). Any 
weight going into the upper level at the — a + l)'th 
crossing goes straight on to the final energy E a , since we 
neglect further deflection at the higher crossings. The 
final (S — 0) energies are E a w n(n-2) _^ ^ with 
a = 0,1, , ...N/2. From this picture, the final energy 
is found to be 

PiE N/2 +(l-p 1 )p 2 E N/2 _ 1 + (l-p 1 )(l-p 2 )p 3 E N/2 _ 2 +.... 
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For moderate values of N for which exact numerical evo- 
lution is feasible up to large times, we find this relation- 
ship to work very well at small K (e.g., Fig. [3^)- 

There are two situations, corresponding to distinct 
physical pictures of the excitation process, where Eq. 
(fTj) reduces to compact forms. First, when r is large 
enough that one can use the Landau-Zener formula, but 
small enough that the pi = e~ aiT are close to unity, the 
excitation at the first two crossings deplete the weight, 
since pi,2 ~ 1. Thus only the top two final levels (E N / 2 
and E N / 2 -i) contribute to the final energy. Expanding 
Pi = e~ aiT to linear order, one gets for this intermediate 
regime 
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Here A£i ns t is the residual energy for the instantaneous 
case (r = 0), namely ( V*=o I H 5= ° |&=o ) - E s =°. We 



have used AEi nst . ~ E N / 2 — Eq ~ N 2 /A. (In an instan- 
taneous quench, the weight would all go to the highest 
final level.) One could also extend the above analysis 
by including the third crossing and expanding up to t 2 . 
Obviously, including further levels and expanding the ex- 
ponentials up to higher orders in r, one eventually gets 
back the full expression (p}. Note that, although we have 
obtained a linear intermediate behavior in r, this is a non- 
perturpative result since it is is based on Landau-Zener 
probabilities. 

Second, at very large r the pi are small, so that 
Pi < ^P2 < ^--- < ^-Pn/2, because of the bosonic factors i(N — 
z + l). Neglecting Pi>i, one obtains 
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FIG. 3: (Color online.) (a) (N,K,S ) = (10,0.1,300). Thick 
line: exact numerical values. The three dashed lines are the 
analytical results for small, intermediate and large r (Eqs. [S] 
[2j [3j) . Eq. Q]is dotted line interpolating between [2l [3l but is 
barely visible because the exact curve almost coincides. (The 
exact curve has some additional oscillations.) (b) Quadratic 
behavior, shown here for cases we do not treat analytically, 
large K (upper) and 80 < NU (lower), (c) Coefficient 77 fol- 
lows Eq. [5] (horizontal line) when Sq^NU. 



Only the lowest and highest energy levels contribute in 
this case. 

Fig. [3^, shows an example where the behaviors of 
Eqs. (P),©,© can be seen in exact numerical calcu- 
lations. The small oscillations at large r on top of 
Eqs. CE]),© are interference effects, discussed later. 



III. FAST QUENCHES 

We now consider small t, i.e., almost instantaneous 
quenches. The main observation is that the residual en- 
ergy has the dependence 
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The quadratic behavior is very robust, and is present 
for all values of N, Sq, K we have checked (Fig. [3^,b). 
In Fig. [3^, we can see the quadratic behavior (straight 
line in log-log plot) for a K<tlNU<^8o case, which is the 
parameter region we analyze. Fig. [3Jd shows quadratic 
behaviors for cases outside this parameter regime. 

To explain this behavior, we start with N = 2 bosons. 
Writing the wavefunction as \ip{t)) = Ci(t)\2, 0) + 
c 2 (t) \1,1) + C3(£) |0, 2), one can consider equations of 
motion for Cj, e.g., c'i(t) = — i [U — S(t)] ci(i) + i-^c 2 (t), 
and solve for small orders in time: 



Ci{t) = Ci(0) + di(0)t + ±Q(0)i 2 + 0(t 3 ) . 



FIG. 4: (Color online.) (a) Quantum interference at large K. 
(b) Interferences resulting from using | AT, ) as initial state. 



Since Ci(0) is entirely imaginary, the linear terms in 
t are imaginary. Linear-i terms will thus cancel out 
from observable quantities like |cj(£)| 2 and C\{t)*c\(t) + 
C\{t)c\{t)* , which appear in the expression for the final 
energy: U (|ci| 2 + |c 3 | 2 ) - ^ [cjc 2 + c* 2 c 3 + h.c.]. The 

energy at r thus has a constant and a t 2 term, but no 
0(t) term. 

For N = 2, we can express the coefficient ?/ in Eq. (|4| 
analytically in terms of the initial Ci(t — 0), and also 
write analytic expressions for Ci(0), as solutions of cubic 
polynomials. Unfortunately, these expressions are too 
cumbersome to be useful. For So ^> U, K, one can calcu- 
late perturbatively, yielding r\ ~ K 2 So/ 24. 

For N > 2, the same argument holds for a leading 
t 2 correction. For very large 6q and small enough r, one 
can use the approximation that only the most imbalanced 
configurations are excited, and restrict to the subspace 
\N,0), |JV — 1, 1>, |JV-2,2). The calculation is then 
similar to the N = 2 case. One obtains 
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at leading order in S . This also contains the N = 2 
result. Fig. [3^,c show that this expression works well for 
S > NU > K. 

While this compact result is valid only for 
5o^>NU^>K, the level crossing structure plays no 
role in this analysis, so in principle it can be applied also 
to other parameter regimes; however the expressions for 
77 are too complicated to be useful in such cases. 

Fig. [3^ also displays the entire instantaneous-to- 
adiabatic crossover, which includes our slow-quench and 
fast-quench results but also a range of r (~0.02 to ^20) 
for which we do not have simple descriptions. 



IV. QUANTUM INTERFERENCES 

Our treatment of multi-crossing configurations (Fig. [2]) 
utilized unit excitation probabilities at the higher cross- 
ings, which is valid at small K. However, if there is 
substantial splitting at each crossing, the configurations 
offer rich possibilities for quantum interference of differ- 
ent paths. We have already seen interference signatures 
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FIG. 5: (Color online.) Fate of self-trapping, (a-d) Num- 
ber imbalance dynamics during and after quench, for various 
quench times. Note different scale on top panel a. (e) Long- 
time average of relative number imbalance after quench is 
over at t = t. Averaging is performed between t = r and 
t — t up = 500000. (f) Quadratic deviation of {z)t> T from 
instantaneous case, for small r. 



in the mild oscillatory behavior in Fig. [3^ for (N,K) = 
(10,0.1) in the near-adiabatic regime. In Fig.[3Ji we show 
more pronounced interference effects at large K. 

Interference between paths also becomes prominent 
when the initial state is not the ground state but has 
weights in more than one eigenstate, e.g., if one starts 
with the state |iV, 0) (Fig. gp). This is not exactly the 
ground state at 5 = — S as long as So is finite. It could 
however easily be an experimentally relevant initial state. 



V. SELF-TRAPPING WITH FINITE-RATE 
QUENCHES 

We now consider the effect of finite quench times on the 
quantum self-trapping phenomenon^ - — Self-trapping 
involves dynamics starting with a biased state, which is 
the case for an instantanous r — quench from large neg- 
ative S = —So- This is the type of dynamics analyzed in 
Ref . [l3l. The relative number imbalance z = (m — n 2 ) /N 
can oscillate around a nonzero value of z if NU/K and 
the starting z are large enough. Since the state ipt=r af- 
ter a finite-r quench is closer to the unbiased [z = 0) 
ground state at S = 0, self-trapping is clearly weakened 
at larger r. 

We can characterize self-trapping through the long- 
time average of z = [n\ — n 2 ) /N at times t > r. A 
nonzero {z) r>T indicates self-trapping. The definition is 
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somewhat subtle because quantum self-trapping involves 
nonzero (z) t y T up to some very large time but not re- 
ally infinite times. (This very large timescale is due to 
the tiny splitting energies of the higher level crossings in 
Fig. [2^,J^1S) The long-time average therefore depends 
on the upper time limit t up up to which z(t) is averaged, 
and will always be zero for large enough t up . This is 
in contrast to the classical self-trapping observed in the 
two-site discrete nonlinear Schrodinger equation, where 
the 'self-trapped' states are truly stationary states^ 

In Fig. [5^-d we demonstrate the weakening of the 
self-trapping effect at larger r, through the behavior of 
z(t). The logarithmic scale highlights the fact that there 
are several imporant time scales in the dynamics. The 
intermediate-r cases have more complicated wavcfunc- 
tions at the end of the quench, as analyzed in Section [TT1 
Hence the dynamics is most interesting for intermediate 
t. This can be seen for example through a larger number 
of peaks in the Fourier transform of the intermediate- 
t time evolution data. In comparison, the r < 1 and 
t 100 cases show fewer high-frequency components in 
the dynamics. 

Fig.JSJs shows the dependence of (z) r <t<t up on r. The 
averaging displayed in Fig. [5^,f was performed up to 
t up = 5 x 10 s ; the z(t) behaviors in Fig. [Sh.-d suggests 
this is a reasonable value for capturing the self-trapping 
phenomenon. The behavior of (z) is similar to the r- 
dependence of the residual energy AE (Fig. [Ip); how- 
erver the ambiguity in the averaging procedure makes it 
difficult to perform a systematic analysis similar to what 
we have presented for AE in previous sections. The sim- 
ilarity of Fig. QJ and Fig. [5]3 is not surprising because 
both the z(t > r) dynamics and the final energy are de- 
termined by the state at the end of the quench, i.e., the 
weight of excited states in ipt=r- At small r there is also 
a quadratic deviation of (z) t>T from the instantaneous 
case (Fig. [Sf), as there is for AE. 

VI. SUMMARY & OPEN ISSUES 

This work analyzes finite-rate quenches most relevant 
to the self-trapping phenomenon, namely, from S = —5q 



to 5 = with parameters -?T<C-/V£/<Ci5o- While quenches 
of the tilt <5 are natural for the dimer, it is relatively new 
in the non-equilibrium literature because bias quenches 
do not appear naturally in the many-site case. We pro- 
vide concise analytical results for the residual energy in 
parametrically different regimes of small, intermediate, 
and large r (Eqs. d [3]). The small r result (JSJ) is 
perturbative in r, while Eqs. [3J [3] are obtained from the 
Landau-Zener formula and are thus non-perturbative, de- 
spite the linear dependence © at intermediate r. 

The importance of Bose-Hubbard dimer dynamics 
reaches far beyond the cold-atom context, as the model 
appears in diverse areas of physics. For example, it is 
equivalent to a large-spin Hamiltonian: H = —KJ X + 
UJz+5J z . This describes single-molecule magnets. How- 
ever for molecular magnet experiments the anisotropy 
(J^) term usually has negative coefficient, and effects of 
host lattice vibrations or nuclear spins cannot always be 
neglected. 

The residual energy is in principle experimentally mea- 
surable in cold-atom realizations through time-of-flight 
measurements that provide energy information from the 
cloud expansion rate (e.g., Ref. \2w. For molecular mag- 
net experiments energy is difficult to measure, but mag- 
netization dynamics (analogous to our z(t) dynamics) is 
commonly reported. 

Our work raises a number of open issues. There are 
a number of parameter regimes other than ours which 
might be of interest. Examples are other values of ini- 
tial and final 5, or of K/U. Refs. [HH have consid- 
ered sweeps of 5 from negative to positive infinity, for 
parameter regions and initial states quite different from 
ours. Clearly, we are only seeing the beginning stages of 
an emerging unified dynamical picture. Another intrigu- 
ing issue is the connection to the mean-field description 
via the discrete Gross-Pitaevskii equation. Ref. [28| has 
worked out some Landau-Zener issues for the mean-field 
dimer, for U < and (-co — > +oo) sweeps. Numeri- 
cal explorations for (— So —> 0) quenches have shown us 
behaviors similar to what we have presented for the full 
quantum case. Finally, the present results may need to be 
adapted for specific realizations, once such experiments 
are designed. 
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